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Abstract 

Correlated stability conjecture (CSC) proposed by Gubser and Mitra [1,2] linked the 
thermodynamic and classical (in) stabilities of black branes. The classical instabilities, 
whenever occurring, were conjectured to arise as Gregory-Laflamme (GL) instabilities 
of translationally invariant horizons. In [3] it was shown that the thermodynamic insta- 
bilities, specifically the negative specific heat, indeed result in the instabilities in the hy- 
drodynamic spectrum of holographically dual plasma excitations. A counter-example 
of CSC was presented in the context of black branes with scalar hair undergoing a 
second-order phase transition [4]. In this paper we discuss a related counter-example 
of CSC conjecture, where a thermodynamically stable translationally invariant horizon 
has a genuine tachyonic instability. We study the spectrum of quasinormal excitations 
of a black brane undergoing a continuous phase transition, and explicitly identify the 
instability. We compute the critical exponents of the critical momenta and the fre- 
quency of the unstable fluctuations and identify the dynamical critical exponent of the 
model. 
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1 Introduction 

Black holes with translationally invariant horizons (also referred to as black branes) 
are ubiquitous in string theory [5]. They are particularly important in the context of 
holographic gauge theory/string theory correspondence [6], where they describe the 
thermal equilibrium states of the dual gauge theory plasma. In case of d = p + 1 
space-time dimensional non-conformal gauge theory plasma, the dual black branes 
of the (d + l)-dimensional supergravity are rather complicated, and have multiple 
scalar hair, see for example [7, 8] . Being dual to a finite temperature gauge theory 
plasma, these black branes have nonvanishing Hawking temperature and necessarily 
are non-supersymmetric. Thus, their classical stability is not assured. The linear 
stability analysis of black branes is complicated by mixing the scalar hair and the metric 
fluctuations. A much simpler exercise is to compute the thermodynamic potentials of 
the black branes, i.e., the free energy density, the energy density, etc., and analyze a 
thermodynamic stability of the resulting potentials. It was a very welcome conjecture 
by Gubser and Mitra [1,2] ■ - Correlated Stability Conjecture — that identified the 
thermodynamic stabilities or instabilities of translationally invariant horizons with the 
classical stabilities or instabilities correspondingly. Moreover, the classical instabilities, 
whenever occurring, were conjectured to appear as Gregory-Laflamme instabilities [9, 
10]. Recall that GL instabilities of the translationally invariant horizons are linearized 
fluctuations <5$ oc e~ l0Jt+lk ' x with a dispersion relation (to, q) such that 



Im(tu) 

where 



> 0, (1.1) 

q<q c 



W \k\ ,„ ^ 

and T is the temperature. Note that we always assume that Im(q) = 0, i.e., fluctuation 
amplitudes do not grow exponentially along the translational directions of a horizon. 
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A standard claim in classical thermodynamics 1 is that a system is thermodynam- 
ically stable if the Hessian mfg A of the energy density £ = £(s,Qa) with respect to 
the entropy density s and charges Qa = {Qi, • • -Qn}, 



a 2 e d 2 s 



ww£ _ ( ds 2 8s8Qb \ (-1 q\ 

M s,Qa= I d 2 S d 2 £ I ' sL " 

\ 8Qa9s 8Qa9Qb 

does not have negative eigenvalues. In the simplest case n = 0, i.e., no conserved 
charges, the thermodynamic stability implies that 

d 2 £ T 

that is the specific heat c v is positive. In the context of gauge theory/string theory 
correspondence black holes with translationary invariant horizons in asymptotically 
anti-de-Sitter space-time are dual (equivalent) to equilibrium thermal states of certain 
strongly coupled systems. Thus, the above thermodynamic stability criteria should be 
directly applicable to black branes as well. CSC asserts that it is only when the Hessian 
(1.3) for a given black brane geometry is positive, the spectrum of on-shell excitations 
in this background geometry is free from tachyons. 

CSC is trivially true in a holographic context, whenever the specific heat of the 
corresponding black brane is negative [3] . Indeed, in the absence of conserved charges, 
the speed of hydrodynamic (sound channel) modes in the dual plasma is 

<£ = f, (1-5) 

c v 

where 5, S are the entropy and the energy densities, and c v 

a> 

is the specific heat. Thus whenever c v < 0, the sound modes in plasma, and corre- 
spondingly the dual quasinormal modes in the gravitational background, are unstable 2 . 

The situation is more complicated in holographic examples in the presence of con- 
served charges. Consider a strongly coupled M = 4 SU(N) supersymmetric Yang-Mills 
theory in the planar limit at finite temperature T and a chemical potential \x for a sin- 
gle global U(l) C 5*0(6) R-symmetry. The dual gravitational geometry is that of the 



1 Assuming that the temperature is positive. 

2 Such instability occurs in N = 2* holographic plasma [7, 11]. 
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Reissner- Nordstrom asymptotic AdS$ black brane [12-14] . For any temperature larger 
than the critical one, T c = fj,y/2/n, there are two phases of the plasma. Although 
one of the phases has a negative specific heat, the sound modes in plasma are always 
stable [15]: 

where 

= v ^+ 1.8 

(J, y/K 

so that given T > T c there are two possible values of k distinguishing two different 
phases: one with k < 2 (a thermodynamically stable phase) and one with k > 2 (a 
thermodynamically unstable phase). There is no contradiction here with the argument 
presented in [3] - the speed of sound in the presence of conserved charges takes a 
more complicated form than in (1.5). In case of a charged plasma we have 

where P is the pressure. Remarkably, there is no contradiction with CSC conjecture 
as well! Indeed, it was shown in [15] that a two-point retarded correlation function of 
the charge density fluctuations has a pole (and thus the SYM plasma has a physical 
excitation with this dispersion relation) at 

to = A i(k - 2) q 2 + £>(q 4 ) , .4 = 0.33339(2), (1.10) 

so that a phase with k > 2 is classically unstable, at least with respect to fluctuations 
with sufficiently small q. We emphasize that even though the gravitational instability 
is in the quasinormal sound channel mode, it is not a hydrodynamic sound, which is 
characterized by the dispersion relation to oc q, see (1.7). 

So far, the examples we discussed provide support for CSC conjecture. It is believed 
however, that CSC conjecture is false [4] 3 . The idea advocated in [4] is to consider 
gravitational backgrounds with translationally invariant horizon that are dual to gauge 
theory plasmas undergoing a continuous phase transition. In the vicinity of the phase 
transition the condensate does not noticeably modify the thermodynamics, and thus 
should not affect the thermodynamic stability of the system. On the other hand, the 



3 See also [16]. 
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phase of the system with the higher free energy is expected to be classically unstable. 
The condensation of the tachyon should bring the system to the equilibrium phase with 
the lowest free energy. The tachyon of the unstable plasma phase should appear in a 
dual gravitational description as a tachyonic quasinormal mode, i.e., the GL instability 
(1.1). 

In this paper we discuss an explicit example of the physical scenario suggested in [4] 
from the perspective of the quasinormal spectrum. We identify the GL tachyon and 
compute the critical exponents for the critical momentum and the frequency of the cor- 
responding fluctuations. Our model is the 'exotic hairy black hole' introduced in [17]. 
The dispersion relation for the sound waves in this model was studied in [18]. Much 
like in case of the R-charged M = 4 SYM plasma, there are no instabilities in hydro- 
dynamic modes. We review the thermodynamics and the hydrodynamics of the model 
in section 2. In section 3 we identify non-hydrodynamic quasinormal excitations in the 
sound channel of the gravitational background with the dispersion features identical 
to those of the GL instabilities (1.1). The GL tachyons exist in the symmetry broken 
phase only, and disappear from the spectrum at the second-order phase transition. We 
compute critical exponents associated with vanishing of q c and ro c = to(q = 0) at the 
transition. Finally, we comment on the dynamical universality class of the model. 

2 Thermodynamics and hydrodynamics of exotic black branes 

The model considered here is not a string theory derived example of gauge/gravity cor- 
respondence — rather, it should be viewed as a phenomenological model of holography. 
Of course, whether or not the model can be embedded into the full string theory is 
irrelevant in so far as one is interested in CSC. The holography here is simply a useful 
tool to think physically about mathematical problem of (in) stabilities of translationally 
invariant horizons. 

Following [19], consider a relativistic conformal field theory in 2+1 dimensions, 
deformed by a relevant operator O r : 

He ft — > H — He ft + \O r . (2.1) 

Such a deformation softly breaks the scale invariance and induces the renormalization 
group flow. We further assume that the deformed theory H has an irrelevant operator 
Oi that mixes along the RG flow with O r . The explicit holographic model realizing 
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Figure 1: (Colour online) The free energy densities Qq of the symmetric phase (red 
curve, left plot) and of the symmetry broken phase (purple curve, left plot) as 
a function of the reduced temperature |r in the gauge theory plasma dual to the 
holographic RG flow in [17]. The right plot represents the square of (Oi) (which we 
use as an order parameter for the transition) as a function of the reduced temperature. 
The dashed green line is a linear fit to (Oi) 2 . 



this scenario was discussed in [17]: 

S, =S CFT + S r + Si= ^ J rfx 4 v^7 [£ C FT + C r + &] , (2.2) 

Ccft = R + 6, L >- = -\ W) 2 + ^ > £i = ~\ ( V ^) 2 " 2X 2 - #Y , (2-3) 

where we split the action into (a holographic dual to) a CFT part Soft', its deformation 
by a relevant operator O r \ and a sector Si involving an irrelevant operator Oi along 
with its mixing with O r under the RG dynamics. The four dimensional gravitational 
constant k is related to the central charge c of the UV fixed point as 

192 . . 

c = — . 2.4 

K 

In our case the scaling dimension of O r is 2 and the scaling dimension of Oi is 4. In 
order to have asymptotically AdS^ solutions, we assume that only the normalizable 
mode of Oi is nonzero near the boundary. Finally, we assume that g < in order to 
holographically induce the critical behavior [19]. 

Effective action (2.2) has a Z 2 x Z 2 discrete symmetry that acts as a parity trans- 
formation on the scalar fields <fi and \. The discrete symmetry <fi — > —cf) is softly broken 
by a relevant deformation of the Ad CFT; while the x ~ >* — X symmetry is broken 
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Figure 2: (Colour online) The speed of sound c s and the ratio of bulk-to-shear viscosities 
in the gauge theory plasma dual to the holographic RG flow in [17] as a function of 
the reduced temperature -Sr. 



spontaneously. Spontaneous breaking of the latter symmetry is caused by the develop- 
ment of the condensate for the irrelevant operator 0\ . The unusual part of this phase 
transition is that the symmetry broken phase occurs at high temperatures (rather than 
at low temperatures) and that the broken phase has a higher free energy density than 
the unbroken phase with (Oi) = 0. The free energy densities of the Z2-symmetric 
phase ((Oi) = 0) and the symmetry broken phase ((Oi) ^ 0) are presented in Figure 
1. 

The hydrodynamic sound channel quasinormal modes of both phases were studied 
in details in [18]. There is a jump discontinuity at the transition in the speed of sound; 
with the speed of sound in the symmetry broken phase being ~ 1% higher than in 
the symmetric phase. Both the shear and the bulk viscosities in the two phases are 
positive. Figure 2 collects the results for the speed of sound and the ratio of bulk-to- 
shear viscosities in the exotic plasma. In the symmetry broken phase the bulk viscosity 
diverges in the vicinity of the phase transition [18] 4 

c 



n 



oc \{Oi 



1-2 



oc (T-T c 



(2.5) 



broken 



Note that since c 2 s > 0, the specific heat (see (1.5)) of both the symmetric and the 
symmetry broken phases is positive — the two phases are thermodynamically stable 5 . 

4 The shear viscosity is continuous and finite at the transition: rj/s = 1/Att [20-22]. 



We also verified that the susceptibility |^ 

s— constant 

the relevant operator O r , see [17] for more details. 



is positive — A is the scale introduced by 
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Figure 3: (Colour online) Dispersion relation of quasinormal modes in the symme- 
try broken (classically unstable) phase of the exotic black branes. The red line (left 
plot) represents the modes at the threshold of GL instability: (to, q) = (0,q c ). The 
green/blue dots (left plot) are the stable/unstable quasinormal modes with ro = — O.li 
and to = O.li correspondingly. The right plot presents the modes at the threshold of 
instability in the vicinity of the critical point. The dashed green line (right plot) is the 
best quadratic in (l — -p 1 ) fit to data. 

The dispersion relation for the sound mode takes the form 

tv = ± Cs q- t - (i + £) q 2 + ^(q 3 ), (2-6) 

implying that the sound modes are stable as well. 

Since the free energy of the symmetry broken phase is bigger, this phase is not 
thermodynamically preferable and one expects a classical instability driving the system 
to a true (symmetric) ground state. Such instability is expected to be of GL type, see 
(1.1), and thus is not expected to show up in the hydrodynamic limit. In the next 
section we explicitly identify the Gregory-Laflamme instability, and thus prove that 
CSC is false. 



3 GL tachyons of exotic black branes 

Physical excitations in a gauge theory plasma are dual to quasinormal modes of a 
corresponding black brane solution [23]. For a translationally invariant horizon, a 
wave-function of a generic background field fluctuation 5$ takes form 

5$ = F{ r )e~ MU , (3.1) 
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Figure 4: (Colour online) Frequency of the homogeneous and isotropic unstable quasi- 
normal modes (tn, q) = (tn e , 0) in the symmetry broken (classically unstable) phase of 
the exotic black branes. The dashed green line (right plot) is the best quadratic in 
(l — ^f) fit to data in the vicinity of the critical point. 

where J 7 is a radial wave function. Fluctuations of different fields will typically mix with 
each other. It is convenient to decompose the fluctuations in irreducible representations 
with respect to rotations about the fc-axis. In the case of the exotic black branes, the 
GL instabilities are expected to arise as fluctuations in the symmetry-breaking (scalar) 
condensate. Thus, they appear in the scalar quasinormal mode channel. It is exactly 
the same set of fluctuations which in the hydrodynamic limit (to — > 0, q — > and 
— — > constant) describe the propagation of the sound waves in plasma. The technique 
for computing the scalar channel quasinormal modes in non-conformal holographic 
models was developed in [24]. It is straightforward to generalize the method of [24] to 
the computation of the scalar channel quasinormal modes of the gauge theory plasma 
dual to the holographic RG flow [17]. The details of the latter analysis will appear 
elsewhere [25] and here we report only the results, focusing on the GL mode 6 . 

The red solid line (left plot) in Figure 3 presents the dispersion relation for the 
quasinormal modes in the symmetry broken (classically unstable) phase at the thresh- 
old of instability: 

(t*,q) = (0,q c ), (3.2) 

threshold 

as a function of the reduced temperature ^. The modes with q 2 > q 2 attenuate, i.e., 
they have Im(tn) < 0, while those with q 2 < q 2 represent a genuine GL instability - 



'Sound waves were discussed in [18] and reviewed in the previous section. 
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they have Im(tt)) > 0. The green dots on the left plot have to = —O.li and the blue 
dots have to = O.li. The right plot on Figure 3 represents the dispersion relation of the 
quasinormal modes in the symmetry broken phase at the threshold of instability in the 
vicinity of the critical point. The green dashed line is the best quadratic in t = (1 — ^) 
fit to the data: 



q c 



= 3.(4) x 10~ 5 + 24.(3) t - 369.(1) t 2 + 0{f) . (3.3) 

fit 



The data suggests that 

q c 2 oc (T — T c ) , (3.4) 

in the vicinity of the critical point. 

The left plot in Figure 4 represents the unstable homogeneous and isotropic quasi- 
normal modes, i.e., GL instabilities with q = 0. The right plot shows the frequency 
dependence of these modes in the vicinity of the critical point. The green dashed line 
is the best quadratic in t fit to the data: 



ito, 



= -9.(8) x 1(T 6 - 7.(2) t + 102.(9) t 2 + 0{f) . (3.5) 

fit 



The data suggests that 

tx» c oc i(T - T c ) , (3.6) 

in the vicinity of the critical point. 

Notice that since Yo c oc i(\ 2 in the vicinity of the critical point, it is natural to identify 
the dynamical critical exponent of the model with z — 2. The symmetries of our model 
identify it as 'model A' according to classification of Hohenberg and Halperin (HH) [26] , 
assuming a natural extension of the classification to classically unstable phases. Mean- 
field models (with vanishing anomalous static critical exponent) in the HH universality 
class A are predicted to have the dynamical critical exponent z mo( i e i-A = 2, as the one 
we obtained 7 . 

In this section we explicitly identified the Gregory-Laflamme instabilities in the 
symmetry broken phase of the exotic black branes. Since this phase is thermody- 
namically stable, our model presents a counter-example to the correlated stability 
conjecture. 



7 The second-order phase transition in exotic black branes is of the mean-field type [25] . 
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